ON THE WEAK-COUPLING LIMIT 
AND COMPLETE POSITIVITY 
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We consider two non-interacting systems embedded in a heat bath. 
If they remain dynamically independent, physical inconsistencies are 
avoided only if the single-system reduced dynamics is completely posi- 
tive also beyond the weak-coupling limit. 
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1. Introduction 

In a variety of different contexts, ranging from quantum optics to the foundations 
of quantum mechanics, the dissipative and irreversible time-evolutions of open quantum 
systems in weak interaction with suitable, large environments are commonly described by 
the so-called quantum dynamical semigroups [1-3]. These consist of linear maps At, t > 0, 
that act on the physical states (density matrices) p s of the open systems, S, and satisfy 
an evolution equation of Kossakowski-Lindlab form [4-10]. 

The maps A t are linear and completely positive, that is, if S is coupled to an arbitrary 
iV-level system Sn, the maps A t <g> id^v preserve the positivity of the states p s+SN of the 
compound system S + Sn for all N. 

Complete positivity is an algebraic property [10] whose physical implications are better 
understood in the negative [11]: if A t is not completely positive, then for some N an 
entangled initial state of S + Sn surely exists which develops negative eigenvalues under 
the action of A t <8> idjv- On the other hand, if the initial state of S + Sn is not entangled, 
i.e. p s+Spf = p s p SN , it never develops negative eigenvalues under the action of At <g>idjv. 

Physical consistency demands that the eigenvalues of At [p s ] must be positive in order 
to be interpretable as probabilities. It is logically necessary, but physically less compulsory, 
that the same should be true of all A t <g> idN[p s+SN ]; this latter request would only be 
guaranteed by the complete positivity of At. Yet, the system Sn is totally arbitrary 
and unchanging, only statistical correlations with S being allowed. Although such an 
occurrence is always possible, it is not always accepted [12] as a justification [11] of complete 
positivity as a necessary property of reduced dynamics. 

It is sometimes argued [12-14] that complete positivity is the consequence of two auxil- 
iary technical simplifications that are essential in the standard [1-3] derivation of quantum 
dynamical semigroups from the closed dynamics of the system S plus its environment. It 
is in fact assumed a) that the initial state of S be uncorrelated to that of the environment, 
and b) that a Markov approximation is possible on rescaled times r = X 2 t, where A << 1 is 
the strength of the system-environment interaction [2-9] (the weak-coupling limit). How- 
ever, requests a) ad b) are not always physically plausible [12-21]; in particular, it might be 
necessary to examine the subsystem dynamics on times of the order of \H, hence beyond 
the weak-coupling limit [15, 17-20] and, in such instances, dynamics not of completely 
positive type may appear. 

The standard derivation of the Redfield-Bloch equations [22], commonly used to de- 
scribe the reduced dynamics of open two-level systems in chemical physics, fails to produce 
even positive dynamical maps A t [4]. While the danger is fully acknowledged [16], it is not 
accepted that one should end up with a completely positive time-evolution [4]. Rather, it is 
argued that only those states whose positivity is preserved should be physically admissible 
or that a slippage in the initial conditions is needed in order to avoid inconsistencies [13-17]. 
However, from the previous considerations, it is clear that accommodating the problem 
of positivity does not properly address the issue of complete positivity which is strictly 
related to quantum entanglement. 

The relevance and role of complete positivity is most clearly seen in the phenomenol- 
ogy of neutral K-mesons as open quantum systems in interaction with a gravitational 
background. Geometrical fluctuations at Planck's scale act as a source of dissipation and 
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decoherence [23-29]. Each single K-meson is thus assumed to evolve according to a semi- 
group of positivity-preserving, entropy increasing phenomenological linear maps A t . It 
turns out [25-29] that these maps must also be complete positive. Otherwise, physical 
inconsistencies would plague the resulting phenomenology of couples of K-mesons evolv- 
ing in time according to the factorized dynamical maps A t <S> A t . More important, these 
dissipative phenomenological models can actually be put to test in experiments performed 
at the so-called (^-factories [26,29]. 

In the following, we show that if the environment is such that there is no induced 
interaction between two otherwise non-interacting systems embedded in it, then the two- 
system reduced dynamics is in factorized form. Moreover this will be true at the 2-nd 
and 4-th order in the system-environment coupling constant A. It then follows that if, 
beyond the weak-coupling limit, the reduced dynamics is not completely positive, then 
either the environment establishes a dynamical dependence between the two subsystems 
or the approximations leading to the reduced dynamics are not physically consistent. 



2. Complete Positivity. 

Let S be a physical system whose time-evolution is given by a (semi) group of linear 
maps At acting on the states of S represented by density matrices p s . Usually, A t is 
positive and thus preserves the positivity of the eigenvalues of p . ^From an abstract 
point of view [10,11], A t is completely positive if and only if the map A t <g> id;v acting on 
the states p s+SN of the system S coupled with an arbitrary iV-level system Sn, is also 
positive for all possible N. 

If A t is only positive, troubles are expected when S is coupled to a generic iV-level 
system and the joint state p s+SN carries correlations between S and Sn- In fact, if S and 
S n are not entangled, then p s + SN = p s <g> p SN (or a convex combinations of factorized 
states), so that 

A t ®id N [p S ®P Sn ] =A t [p s ] 



-.Sn 



(2.1) 

and the positivity of A t <E> id^v automatically follows from the positivity of A t . 

However, let S be a two-level system and Sn = S 2 = S. As a common state p s+s we 
consider the projection 
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(2.2a) 
(2.2b) 



onto a singlet-like state of S + S with eigenvalues and 1. As a linear map on S, let us 



consider the transposition operation T 
but 
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The map T is positive, 



T®id 2 [p A ] 
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(2.3a) 
(2.3b) 
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has eigenvalues ±1/2. Therefore the transposition T is not completely positive, already 
the coupling to a 2-level system failing to be positive [5] . Clearly, the origin of troubles is 
the term (2.2b) which encodes the entanglement between the two systems and is changed 
by T ® i&2 into (2.3b). Instead, the term (2.2a), which represents an uncorrelated density 
matrix, is left unchanged. 

Stinespring's decomposition theorem [11] ensures that the most general completely 
positive linear map A on the states of S is of the form 

A[p S ]=J2 V iP Sv ^ ( 2 - 4 ) 

i 

where Vg are suitable bounded operators. Evidently, A <g> idArfp 5 " 1 " 5 ^] is positive for any 
positive p s+SN . The quantum mechanical time-evolution 

A t [p s ] = e~ itH p s e itH ' (2.5) 

is of the form (2.4) and maps pure states into pure states thus preserving coherence. 
However, if one considers open quantum systems, coherence is usually lost. Hence, A t 
cannot be of the form (2.5); whether it is of the form (2.4) must be decided on physical 
grounds. 

Remark 2.1 In the approach to if-mesons as open quantum systems [23-26], the trans- 
positio n map T is replaced by a phenomenological dynamical map At and, in [24,26], 
T <S> id.2 by At ® At. It can be shown [26,27,29] that At must be completely positive. 
Otherwise, physically realizable initial correlated states of two K-mesons as in (2.2) would 
develop negative eigenvalues. 

3. Reduced Dynamics 

The quantum open system of interest S is assumed to be weakly interacting with a 
large (infinite) environment R, the dynamics of S + R being governed by the Hamiltonian 

H = H S + H R + XH SR , (3.1) 

where Hs and H R are the Hamiltonians of the system S, respectively environment R and 
Hsr is an interaction term with coupling strength A. The total system S + R is closed 
and its states, represented by density matrices p, evolve reversibly according to 

^ = L H [p t ]:=-i[H,pt] . (3.2) 



The environment is supposed to be in equilibrium with respect to H R , thus 
Let R be a heat bath in equilibrium at temperature namely we take 



P R ,H R 



r _ exp(-(3H R ) 
9 -Trexp (-/?#*) • {d - S) 
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The interaction term is chosen to be of the form 



H S R = J2 A a® A ai (3-4) 

a 

where the self-adjoint operators Af and A R refer to the system and environment, respec- 
tively. It is no restrictive to assume (pn(A R ) = for all a. 

Remarks 3.1 

i) The environment must eventually be considered infinite dimensional in order to allow 
for continuous spectra of Hr and avoid recurrences. The states p R are thus not 
confined to density matrices, nor are the expectation values 

<p R (A R ) = Tr R (p R A R ) (3.5) 



always expressible via a trace operation. However, we will stick to the density ma- 
trix notation, the genuine case of infinitely many degrees of freedom being easily 
retrieved [5-8]. 

ii) The dissipative effects perturb the time-evolution of S given by its own Hamiltonian 
Hs and are at least of second order in the coupling A. The weak-coupling limit 
consists in going from the fast-time variable t to the slow-time variable r = A 2 t, 
with A — > 0. The technical procedure is physically justified when the ratio tr/ts 
between the characteristic time tr of the environment and the characteristic time 
of the dissipative effects on S, rg, is small [3]. However, there might be heat bath 
temperatures for which one has to retain higher powers in tr/ts, being forced to go 
beyond the weak-coupling limit [15]. 

In the interaction representation, p t := exp(— tL Q ) [p t ] with L [p] := —i [Hs + Hr , p], the 
evolution equation reads 

!> "' = A e~ tL ° Lh sr e tL ° [pt] • (3.6) 



dt 

For the sake of simplicity, we assume the spectrum of Hs to be discrete and non- 
degenerate with eigenvalues e r and eigenvectors \r). We enumerate the operators \r)(s\ 
by denoting them as Vf , so that H SR = ^Vf <g> V R with V R = £ Tr s (VfA% )A R . 
Further, setting ujj := e r — e s , we get 

H SR (t) = e- tLo [H SR ] = e "" J 1 V i ®z ltHR V R e~ ltHR . (3.7) 

3 

The formal solution of (3.6) is p t = (l + J2n=i A " U t ] ) [Po], with 

u t n) = dt i / dt 2 - - dt n L HsR ( tl )L HsR ( t2 ) ■ ■ ■ L HsR{ t n) , (3.8) 

Jo Jo Jo 



where La[-] '■= —i 



A, 
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In order to extract the system S reduced dynamics, we operate on the states p of 
S + R with the projector P[p] = Ttr(p) ® Pr which decouples the environment degrees of 
freedom. Further, we take as initial state of S + R the state Pq = Pq = p s ® p R with no 
correlation between S and R. It follows that P[po] = po- 

Remark 3.2 Despite the fact that they are the most used [1-3], the choice of P made 
above and the assumption on the initial global state cannot be generically upheld. In 
particular, one cannot always benefit from a factorized initial state. However, this is in 
many instances plausible, as in the case of neutral K-mesons in a gravitational background. 
Indeed, K- mesons produced in strong <p- meson decays are arguably not influenced by ge- 
ometrical fluctuations of gravitational origin . In general, one may be forced to adopt 
different projectors suited to initial states where system S and environment R result cor- 
related by interactions prior t = [13-17, 21]. 

We now elaborate more in detail on the approach of [19]. We assume the environ- 
ment to be a Bose thermal bath described by the equilibrium state (3.3). The projector 
P involves bath expectations with respect to (3.3), then only even correlation functions 
survive. Keeping terms up to A 4 , one eventually finds 



^ = X>PuPP[p t 



+ X 4 



Ptli A) - ptif ) puf ) 



P[pt 



(3.9) 



where ifj: n ^ is the time-derivative of U^ n \ In particular, the 2-nd and 4-th order contribu- 



tions, PU$ 2) and PUf'PUf' read 



r(2) D rr( 2 ) 



P 



(2) 



dtxP L HsR{ t) L HsR<ytl ) 



(3.10a) 



U^PU\] 2) = j dfc / dt 2 [ t2 dt 3 PL HsR(t) L HsR{tl) PL HsR{t2) L HsR{t3) .(3.10b) 
Jo Jo Jo 

Remark 3.3 After standard rearrangement of the integrals in (3.10b), the whole 4- 
th order contribution in (3.9) assumes a typical cumulant expression [19,20]. Thermal 
correlation functions are expected to factorize for large times; in such a case the cumulants 
vanish and allow one to operate a Markov approximation also at 4-th order in A. 

Since p R commutes with Hr, setting pf := Tr^(p t ), the time-evolution equation 
obeyed by the open system S has the form 



dpi 
dt 



= L fls bf]+A^fVf] + A 4 ™. 



(3.11) 



Let Qf k {t) := <Pr (v^^^Vj^j denote the environment two-point correlation func- 
tions. Then, the 2-nd order dissipative contribution in (3.9) explicitly reads 



K 



(2) 



[p s \ = E / dti e 



-iti uji 



vf P s 



v£ +n; k (ti) 



vc 



P s vi 



(3.12) 
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From (3.10b) it follows that, because of the chosen projector P, the 4-th order dissi- 
pative operator involves four-point thermal correlation functions, which in turn are 
linear combinations of two-point ones. After a lengthy calculation one arrives at 
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K'i'\p 



E E 

j,k,i,m p=l 



JO 



t-ti 



-t 2 



dte~ lA 



jkim 



jkim I 



(3.13) 



where t= (ti,t 2 , t 3 ) and A k£m ■ t= (u k + u e + uJm)t\ + {^t + u m )t 2 + ^mh- 

The quantities QS p \t) are products of two point-correlation functions, while the op- 
erators Dj k \ m are essentially double commutators of observables of S. For instance, in the 
case p = 1, one has 



jkim 



3 



V s p s 

r m r 



(3.14) 



Because of the explicit dependence of both and on time, the right hand side 
of (3.11) retains memory effects and does not generate a semigroup. However, a Markov 
approximation can be performed based on the following argument (see Remark 3.3). The 
two-point correlation functions <p r(V r (t)V R ) are expected to factorize for t larger than 
the correlation-time t r of the environment which is much shorter than the typical time for 
the dissipative effects being felt by the subsystem S. Since we assumed lpr(A r ) = 0, 



it follows that ipii{V-{t)V£ 



for t » t . Therefore, for times t » r , the 



time-dependent dissipative operators and can be replaced by time-independent 
dissipative operators and , by extending to infinity each time-integration in (3.12) 
and (3.13). For a rigorous approach to this kind of Markov approximation the reader is 
referred to [4,6]. Explicitly, 



K^[p s ] := Wo+K) [v k S P S ,Vf] + O-(^) [v k S ,p S Vf 

3,k I 
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j,k,£,mp=l 



(3.15a) 
(3.15b) 



where fi± (u/) := dt e~* u * Sl%(t) and 



/•oo />oo poo 



(P) (f) 



(3.16) 



4. Non-interacting Open Quantum Systems 

Let the open system S consist of two non-interacting systems Si and S 2 whose dy- 
namics, disregarding for the moment the presence of the environment R, is governed by the 
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Hamiltonian operators H$ 1 and H$ 2 - Again, we assume the energies e ar of H$ a , a = 1, 2, 
to be discrete and non-degenerate and set u a j := e ar — e as . 



In absence of the system S = Si + S2 would evolve in time according to the 
Hamiltonian Hs = Hs ± <8> id2 + idi <g> Hs 2 - Instead, we suppose £1,2 to interact weakly 
and independently with a thermal bath in the state (3.4). We chose an interaction term 
of the form 

H S R = E(V5 ® id 2) ® V 5 + ( id l ® ^2?) ® ^2?) • (4-1) 

with ^(V^) = for all j, a = 1, 2, and same coupling constant A. 

The analysis of the previous section can be repeated in this new context, the major 
difference being that, inserting (4.1) in (3.10) extra-indices appear identifying the system 
S a , a = 1,2, the various operators refer to. As a consequence, the dissipative operator 
in (3.12) becomes 



a, 6=1 



(4.2) 

where ^j. hk {t) = <fiu(v^j(±t) V^j and, for sake of simplicity, denotes either Vfj ®id 2 
or idi <S> V^. I 11 turn, the 4-th order operator reads 



10 r t r ti r t-t!-t 2 

-(4) _ / / / Af e -iA bk;ce . dm -t n (p) (I\ D (p) 

ute * L aj;bk;ct,dm\ L ) ^ aj;bk;c£;dm ' 



a1 4) = E E 



(4.3) 



j ,k,l,m 



where A bk;c£ . dm ■ t = (u bk + u c t + u dm )ti + (u c i + u dm )t 2 + ^dmh and 

Xt)--=nt j]ci (tl+t2)ntk;dm(t2 + t3) 



aj\bk\ct\dm * 



(1) 

aj;bk;cl;dm 



[P S ] := 



Vbk > Vet 



vLp s 



(4.4a) 
(4.4b) 



The two-point correlation functions ^aj-bk(^) involve either bath operators interacting 
with the same system, a = c, or with different systems, a ^b. We can perform the Markov 
approximation exactly as in the case of just one system S so that 



K (2) [P S ] :=E|^;a>a,) [^Y > V bk] + Kj;bkM [v& , P S V* 

a,b ^ 
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a ,b,c,d p — ^ 
j ,k ,£,m 



where := /~df e-* w *n± ;6fc (f) and 



(4.5a) 



(4.5b) 



^ij;6fc;c£;dm(^W=;^;dm) 



dfx/ df 2 / dt 3 e- <A "^- r nW fcfc;rf;dra (f) . (4.6) 



Using (4.1) and (3.7), the 2-nd order dissipative operator can be written as follows: 



K[ 2) ®id 2 [p a ] + id! ®K^>[p 



( 2 )r^ 



+ Ei^ W K,) 



ajtb 



V a S iP S , Vr S 



bk 



+ n aj;6feKi) 



* bk i \ J y aj 



(4.7a) 
(4.7b) 



The terms K { a\ a = 1, 2 in (4.7a) are dissipative operators of the form (3.15a), involving 
only observables referring to the system S a - If the contribution (4.7b) were absent, the 
right hand side of (4.7a) would generate a factorized time-evolution p s i— > pf = Aj®A 2 [p s ], 
with K\ satisfying 



id. 



(2) 



id 2 



A, 1 



id 2 k' 



(4.8) 



and analogously for A 2 . 

Clearly, the term (4.7b) dynamically couples the two systems Si and S 2 through their 
interaction with the same environment. However, this coupling depends on the strength of 
the thermal correlations between bath operators describing the interaction with different 
subsystems. If there are no correlatoins between them, that is if 



(4.9) 



whenever a 7^ c and t > 0, then the term (4.7b) does not contribute since, with no 
restriction, we can assume the one-point bath correlation functions to vanish. 

The same result follows from (4.9) when we include 4-th order dissipative effects. This 
can be seen as follows. In the expression (4.4b), the right hand side vanishes if b 7^ c, since 
then V^l and belong to different subsystems and thus commute. With b = c, using (4.9) 
one sees that the right hand side of (4.4a) vanishes unless a = c = b = d. Thus, if (4.9) 
holds, only a = b = c = d = l and a = b = c = d = 2 contribute to the sum in (4.5b). 
Then the 4-th order dissipative operator splits as the 2-nd order one, 



K^[p s ] = K{ 4) ® id 2 [p s ] + id! <g> K [ 2 4) [p 



(4.10) 



where 
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D 



(p) 



K-a ^ — ^aj;ak;ai;ami^cik;at,am ' L '„j:<ik:a< uim ' 

j,k,£,m p=l 



a = 1,2 . 



(4.11) 



As a consequence of (4.9), the evolution equation (3.11), after the Markov approximation, 
becomes 



dtp? = 



(^L Hsi + + K^) ® id 2 + id! ® (l Hs2 + + [pf] (4.12) 
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and the dynamical maps generated by it factorize into Aj ® Af , where A \ and A^ are 
generated as in (4.8) with K [ a \ a = 1, 2, replaced by + . 

5. Conclusions 

The request that the reduced dynamics of open quantum systems in interaction with a 
reservoir be completely positive is often rejected as not physically necessary [12-14,16-21]; 
it lacks physical appeal and involves trivial and uncontrollable couplings of the systems 
with generic iV-level systems. 

The issue of complete positivity assumes its full physical significance when dealing 
with the dynamics of correlated systems interacting with an environment. All depends on 
whether the reduced dynamics of the subsystems factorizes as A t <S>A t , thus indicating that 
they evolve independently from each other, remaining dynamically uncorrelated. 

If this is the case, then the complete positivity of the single system reduced dynamics 
A t is unescapable, otherwise the joint reduced dynamics A t ® A t generates unacceptable 
negative probabilities. ^From this point of view, the necessity of complete positivity ap- 
pears to be a dynamical aspect of quantum entanglement. 

On the other hand, there are situations for which the factorization of the reduced 
dynamics is not the case due to the subsystem-environment interaction and to the physics 
of the environment itself, namely the behaviour of its correlation functions. If the two- 
system reduced dynamics does not factorize, it is not compelling that the single-system 
reduced dynamics be completely positive. 

In the literature there is evidence of non-completely positive reduced dynamics beyond 
the weak-coupling limit [17-20], namely taking into account contributions of order A 4 in 
the coupling between subsystem and environment. We have shown that if the reduced 
dynamics of two non-interacting subsystems factorizes at 2-nd order in A, it factorizes 
also at 4-th order. Therefore, absence of complete positivity in the single-system reduced 
dynamics at 4-th order would jeopardize the description of the physical behaviour of two 
of these subsystems at the same order of approximation. 

Finally, we would like to stress that there are physical instances where an experimental 
check of complete positivity seems achievable, as in the case of X-mesons [25-29] ; there, the 
conditions for a factorized reduced dynamics are plausibly fulfilled because of the weakness 
of the effects of the gravitational background. The reduced dynamics of K-mesons as open 
quantum systems must then be completely positive and experiments at (^-factories can 
explicitly clarify this fundamental request. 
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